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Background: In low-risk gestational trophoblastic neoplasia (GTN) patients, a predictive marker for early identification of 
methotrexate (MTX) resistance would be useful. We previously demonstrated that kinetic modelling of human chorionic 
gonadotrophin (hCG) measurements could provide such a marker. Here we validate this approach in a large independent patient 
cohort. 

Methods: Serum hCG measurements of 800 low-risk GTN patients treated with MTX were analysed. The cohort was divided into 
Model and Test data sets. hCG kinetics were described from initial treatment day to day 50 using: '(hCG(time)) = hCGO*exp(- 
/c'^time) + hCGres', where hCGres is the modelled residual production, hCGO is the baseline hCG level, and k is the rate constant. 
HCGres-predictive value was investigated against previously reported predictors of MTX resistance. 

Results: Declining hCG measurements were well fitted by the model. The best discriminator of MTX resistance in the Model data 
set was hCGres, categorised by an optimal cut-off value of >20.44IUI~^: receiver-operating characteristic (ROC) area under the 
curve (AUC) = 0.87; Se = 0.91; Sp = 0.83. The predictive value of hCGres was reproducible using the Test data set: ROC 
AUC = 0.87; Se = 0.88; Sp = 0.86. Multivariate analyses revealed hCGres as a better predictor of MTX resistance (HR = 1.01, 
P<0.0001) and MTX failure-free survival (HR = 13.25, P<0.0001) than other reported predictive factors. 

Conclusion: hCGres, a modelled kinetic parameter calculated after fully dosed three MTX cycles, has a reproducible value for 
identifying patients with MTX resistance. 



Human chorionic gonadotrophin (hCG) is a sensitive serological 
marker in monitoring trophoblastic disease activity in gestational 
trophoblastic disease (GTD; Bagshawe, 1992). Patients diagnosed 
with low- risk gestational trophoblastic neoplasias (GTN), comprising 
^95% of malignant GTD forms (including invasive mole, 
choriocarcinoma, epitheliod trophoblastic and placental site 



trophoblastic tumours) have International Federation of Gynecology 
and Obstetrics (FIGO) 2000 risk scores ranging from 0 to 6 and 
are treated using single-agent chemotherapy (Kohorn, 2001; 
Chalouhi et al, 2009). The 8-day methotrexate (MTX) protocol 
modified by Bagshawe et a/, 1989 is the most commonly used 
regimen in Europe (Savage et al, 2008; Chalouhi et al, 2009; 
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Kerkmeijer et al, 2009). MTX treatment is continued until either 
normalisation of serum hCG concentration and further three cycles 
or resistance to MTX therapy is detected (Seckl et al, 2010). In cases 
of MTX resistance, observed in nearly one-half of low-risk cases 
(Sita-Lumsden et al, 2012), patients are changed to either 
actinomycin D (ACT-D) or to a combination of etoposide, MTX, 
ACT-D, cyclophosphamide and vincristine regimen (EMA-CO) 
regimens (Powles et al, 2007; Patel & Desai, 2010). Tumour 
chemoresistance is frequently defined as an increase or stagnation 
of hCG levels over a 2- to 3 -week period, but no consensus guideline 
has been defined (Foulmann et al, 2006). Ambiguity concerning the 
definition of MTX resistance based on changing hCG levels increases 
the likelihood that patients will receive repetitive ineffective cycles of 
MTX chemotherapy. 

Several previous studies have attempted to develop new 
methods to enable earlier identification of patients that will 
develop MTX resistance (van Trommel et al, 2006; Savage et al, 
2008; Kerkmeijer et al, 2009; Growdon et al, 2009a). Among these, 
a previous study led by the French reference centre for treatment of 
GTN (Centre de Reference des Maladies Trophoblastiques, Lyon, 
France), showed that mathematical modelling of individual 
declining hCG measurements using a population kinetic approach 
is feasible, and that derived parameters were strong early predictors 
of MTX resistance (You et al, 2010). The study concluded that 
validation of the study's assumptions in independent cohorts of 
patients was warranted. 

Charing Cross Hospital's Trophoblastic Disease Unit (London, 
UK) is a leading centre for treatment of GTN (Seckl et al, 2010). 
The objective of the present study was to validate the predictive 
value of the hCG kinetic modelling approach with respect to MTX 
resistance in low-risk GTN, using the Charing Cross Centre's data. 



PATIENTS AND METHODS 



Patient selection and hCG analysis. HCG measurement and 
clinical data from patients diagnosed with low-risk GTN and 
registered between August 1991 and January 2011 were retrieved 



from the Charing Cross Hospital's database. All patients were 
treated with the 8 -day MTX regimen of 50 mg intramuscular MTX 
on days 1, 3, 5 and 7, combined with 15 mg oral folinic acid on 
days 2, 4, 6 and 8, every 2 weeks as first-line treatment for low-risk 
GTN according to the FIGO 2000 guidelines (Bagshawe et al, 1989; 
Kohorn, 2001). Treatment was continued for three cycles after 
normalisation of serum hCG concentration (<5IUl~^). In cases 
of MTX resistance, defined as three static or three rising hCG 
measurements, treatment was changed to second-line treatment 
(Seckl et al, 2010). The static versus changing distinction being 
± 10% of the preceding serum hCG concentration. 

Exclusion criteria included: (i) treatment with <3 cycles of 
MTX, (ii) treatment variants from the conventional MTX regimen 
and (iii) hCG monitoring with < 3 measurement time points. 

Total hCG serum concentrations were determined using an in- 
house competitive radioimmunoassay, which uses a polyclonal 
rabbit anti-serum. The lower limit of detection is 1.6IUl~^ 
(calculated as 2 s.d. above the zero standard) and the coefficient of 
variation (CV) at ^3IUl~^ is estimated at 14% (Harvey et al, 
2010). 

Model and test patient data sets. The 800 low-risk GTN patient 
cohort treated with the 8-day MTX regimen was divided into two: 
(i) a Model data set composed of 418 patients (195 MTX- resistant 
and 223 MTX-sensitive patients) and (ii) a Test data set 
comprising 382 patients (144 MTX-resistant and 238 MTX- 
sensitive patients; Figure 1). 

Procedure used to ascertain the predictive value of modelled 
kinetic parameters. A three-step scheme was used to ascertain the 
predictive value of population kinetic modelling of hCG data, using 
the patient cohort's Model and Test data sets (Figure 1) according 
to the following outline procedure: 

• Estimation of individual kinetic parameters using population 
kinetic modelling of hCG measurements in Model data set 
patients. Receiver- operating characteristic (ROC) curves were 
used to assess the predictive values of the modelled kinetic 
parameters regarding MTX- resistance risk. 



Whole cohort 

(Low-risk GTN treated with ^ 3 MTX cycles and with ^ 3 hCG time points 
between day 0 and day 50) 
800 patients 



1) Estimation of 
model parameters 
and assessment of 

the predictive 
value of modelled 
kinetic parameters 
for MTX resistance 



2) Independent 
validation of the 
model determined 

and of the ^ 
predictive value of 
hCGres for MTX 
resistance 



3) Comparison of 
hCGres predictive 
value with other 
reported predictors - 
using ROC curves, 
logistic regression 
and MTX failure-free 
survival 



MTX sensitivity : 
223 patients (53.3%) 



MTX resistance: 
195 patients (46.7%) 










Model data set: 
418 patients 





MTX sensitivity : 
238 patients (62.3%) 



MTX resistance: 
144 patients (37.7%) 




© 



Test data set: 
382 patients 



® 



Multivariate analysis data set: 
669 patients: 

MTX resistant, n= 267 (39.9%); MTX sensitive, n= 402 (60.1%) 



Exclusion: 
Patients with 
missing 
data for 
multivariate 
analyses* 

131 patients 



Figure 1. Flowchart of analysed patients in the three-time strategy. *WHO-FIGO score (classified between 0 and 6); choriocarcinoma (yes vs no); 
hCG level measured 1 week after treatment start; hCG level in the seventh week. 
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• Challenging the Test data set using the most predictive kinetic 
parameter identified, using the Model set to estimate the test's 
performance and assess reproducibility. 

• Comparison of the relative predictive value of the best kinetic 
predictor identified using the Model data set with that of 
previously reported predictors of MTX resistance, using ROC 
curves and multivariate statistics (logistic regression; Cox- model 
survival test) in a combined data set, called the multivariate 
analysis data set, comprising model and test patients that had all 
required data (Figure 1). 

Step 1: population kinetic modelling of hCG measurements in 
model data set patients. The population pharmacokinetic model 
developed by the French national centre in 2010 was used to fit 
declining hCG measurements from data of patients extracted from 
the Charing Cross Hospital database (You et al, 2010). Details of 
the mono -exponential model have been previously described (You 
et al 2010). The model was further refined and the kinetic 
parameter 'hCGres', related to residual production at the end of 
hCG decline, was derived following fitting of declining hCG values 
using the final model: 

hCGijit) = (hCGO,-^e"^^'* + hCGres,) (1 + ely) 
where hCGy (t) is the jth hCG measurement in patient / at time t; 
hCGres/ is the hCG residual production of patient /; sly represents 
the proportional residual variability at the jth measurement of 
patient /, which is assumed to follow a normal distribution with 
zero mean and variance. Ki is the hCG decline rate constant. 
Parameters K and hCGres were assumed to vary between patients 
according to a normal log distribution; because of the hetero- 
geneity in baseline hCG levels, a Box-Cox distribution was 
assumed for hCGO in the population. Moreover, censored 
observations of hCG titres below the limit of quantification (i.e., 
those measurements recorded as <2IUl~^) were handled using 
the so-called M3 method, which considers the likelihood of 
observations being below the quantifiable limit (Beal, 2001; Ahn 
et al, 2008; Bergstrand and Karlsson, 2009). The kinetic analysis 
used hCG measurements taken between the first day of MTX 
treatment (day 0) and day 50. If the treatment was altered during 
the first 50 treatments days, further hCG titres were not included 
in the analysis. The predictive performance of the final model was 
evaluated by visual predictive check (VPC) method, an approach 
commonly adopted for internal validation in modelling studies 
(Brendel et al, 2006; Dartois et al, 2007). The analysis was 
performed using a non-linear mixed-effects modelling strategy, 
implemented in NONMEM Version 7 (Beal et al, 2009) (full details 
in Appendix online). 

Assessment of predictive value for MTX resistance using ROC 
curve analysis in model data set patients. Continuous ROC curve 
analyses were used to assess the predictive values of modelled K 
and hCGres in the Model data set. Among them, the hCGres 
kinetic parameter, which offered the highest discriminative value, 
based on the highest value of ROC area under the curve (AUG) was 
selected for further analyses. The optimal hCGres cut-off value for 
discriminating resistant vs sensitive patients was determined by the 
intersection of minimal overlap of the true-positive and true- 
negative populations for the Model data set, that is, where 
sensitivity (Se) + specificity (Sp) was maximised. 

Step 2: validation of 'hCGres' predictive ability using Test data 
set. The final model and optimal cut-off value derived using the 
Model data set was applied, to independently estimate individual 
values of hCGres in the Test data set patients. As previously, 
declining hCG values observed during the first 50 MTX treatment 
days were fitted using the Maximum A posteriori algorithm in 
NONMEMTM software. The theoretical MTX resistance risk of 



the Test data set patients was defined in every subject based on the 
Model set's defined cut-off value. The accuracy of prediction was 
then assessed using discontinuous ROC curve analysis. 

Step 3: assessment of predictive value of hCGres against 
previously reported MTX resistance predictors. Other previously 
reported predictive factors for chemoresistance considered in the 
comparative analysis were: WHO-FIGO score (classified between 
0 and 6); choriocarcinoma (yes vs no; You et al, 2010); hCG level 
measured 1 week after treatment start, categorised according to 
Growdon et al (2009b) cut-offvalue > 2000 lU 1 " ^ (Growdan et al, 
2009b) or hCG level in the seventh week, categorised according to 
different thresholds, ^520.24 vs > 520.24 lUl"^ (van Trommel 
et al, 2006); ^ 500 vs > 500 lU 1 " ^ (Savage et al, 2008); and ^ 737 
vs >737IU1"^ (Kerkmeijer et al, 2009). 

Using the multivariate analysis data set, the relative predictive 
performance of modelled hCGres was compared against the 
previously reported predictors using ROC curve analysis. More- 
over, the independence of the predictive value of hCGres for MTX 
resistance was compared against these predictors, considered as 
continuous or categorised covariates, using multivariate logistic 
regression. A backward elimination procedure was used to select 
the independent predictive factors. 

Predictive value of hCGres regarding MTX failure-free survival 
against previously reported predictors. The data set was also used 
to assess hCGres -predictive ability regarding MTX treatment 
failure. MTX failure was defined as MTX resistance or relapse of 
disease. Relapse was defined as recurrent disease requiring further 
chemotherapy marked by elevated hCG detected following pre- 
viously successful chemotherapy (where hCG had been normalised 
for ^6 weeks). The predictive value of hCGres was assessed using 
univariate (log-rank test) and multivariate (Cox model) survival 
tests. A backward elimination procedure was used to select the 
independent predictive factors. All tests were performed with S-Plus 
(TIBCO Software Inc., Palo Alto, CA, USA) or SPSS (SPSS Inc., 
Chicago, IL, USA) using a two-sided 0.05 a-risk. 

Landmark timepoint analyses. Landmark time point analyses 
with time threshold set at 50 days were performed to ensure the 
validity of results and lack of bias. As a result, ROC curve analyses 
and survival tests were performed with specific data sets, excluding 
patients treated with a second-line treatment before day 50. Results 
are shown in Appendix, (online only) because results were 
consistent with those obtained with all data sets. 



RESULTS 



Patient characteristics in model and test data sets are presented in 
Table 1. 

Population kinetic modelling of hCG kinetics in model data set 
patients. Parameter estimates and their respective s.e. are reported 
in Table 2. Typical values and CV between patients of hCGO, K and 
hCGres were 14 400IUr^ (CV: 152%), 0.169 per day (CV: 306%) 
and 24.7 IU1~^ (CV: 35%), respectively. The goodness-of-fit plots 
show that individual hCG profiles were well fitted during the first 
50 treatment days by the model (Figures 2 A and 3). The VPC 
method showed that the median, 5% and 95% of the observed hCG 
values over time were included within 95% confidence interval 
boundaries of the simulated median, 5% and 95% of hCG, in 
keeping with the good data fitting by the model (Figure 2B). 

Prediction of MTX resistance using derived kinetic parameters 
in model data set patients. The abilities of hCGres and K for 
predicting MTX resistance were compared using ROC curves. On 
the basis of the AUG values, hCGres was a better predictor of MTX 
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Table 1. Summary characteristics of patient cohort 






Model data set 


Test data set 




No. patients 


No. patients 


Age (IQR): years 


31.2 (26.1-36.1) 


30.9 (26.1-35.6) 


Interval between evacuation and 


1.9 (1.3-2.9) 


1.9 (1.4-3.2)^ 


first MTX dosing (IQR): months 






Choriocarcinoma (%) 


19 (4.5) 


18 (4.7) 


WHO-FIGO sore 


0 (%) 


37 (8.9) 


33 (8.6) 


1 (%) 


76 (18.2) 


73 (19.1) 


2 (%) 


95 (22.7) 


82 (21.5) 


3 (%) 


105 (25.1) 


108 (28.3) 


4 (%) 


70 (16.7) 


60 (15.7) 


5 (%) 


28 (6.7) 


18 (4.7) 


6 (%) 


7(1.7) 


8(2.1) 


Number of MTX received: median 


6 (4-7) (3-14) 


6 (4-7) (3-12) 


(IQR) (range) 






Response to MTX 


Sensitive (%) 


223 (53.3) 


238 (62.3) 


Resistant (%) 


195 (46.7) 


144 (37.7) 


Total (n = 800 patients) 


418 


382 


Abbreviations: IQR = interquartile range; MTX = methotrexate; WHO 


-FIGO = WHO-lnter- 


national Federation of Gynecology and Obstetrics. 




^Unknown for one patient. 







Table 2. Typ 

set patients ( 


cal parameters and precisior 
n = 418 patients) 


1^ of estimates for 


Model data 


Parameter 


Population 
value 


Median of 
individual 
predictions 


Interindividual 
variability 


Residual 
variability 


hCGO 

(iui-^)^ 


14 400 (5.7%) 


16 574 


1 52% 


34.4% (1 .4%) 


K (per day) 


0.169 (1.3%) 


0.174 


306% (6.7%) 


ND 


hCGres 

(lui-^) 


24.7 (11%) 


21.97 


35% (4.9%) 


ND 


Abbreviations: hCG = human chorionic gonadotrophin; hCGres = hCG resistant; ND = not 
determined. 

^Precision of parameter estimates were computed from a 100 resampling bootstrap 
analysis. 

'^hCGO individual parameters were assumed to follow a Box-Cox (BOCOX) distribution and 
were implemented as follows: BCOX1 =((EXP(ETA(1)))**SHAPE1— 1)/SHAPE1; hCGO = 
THETA(1)*EXP(BCOX1) with SHAPE1 estimated at -0.207 (s.e.: 7.6%), and ETA(1) a random 
variable following a normal distribution of mean 0 and variance estimated at 4.24 (s.e.: 3.5%). 



resistance than K (AUG for continuous ROC = 0.94, 95% 
CI = 0.91-0.96 vs AUC = 0.69, 95% CI = 0.64-0.74; Figure 4A) 
and was, therefore, selected for further analyses. The optimal cut- 
off value for hCGres for discriminating resistant vs sensitive 
predictor was found to be 20.44 IU1~\ When the model data set 
patients were reassessed by discontinuous ROC analysis after 
classifying them into two predictive groups, as either high MTX- 
resistant risk (> 20.44 lU 1 " ^) or low risk (^20.44 lUl" ^) to 
ascertain hCGres's predictive diagnostic performance at the 
optimal cut-off value, the AUC was 0.87 (95% CI = 0.65-1) with 
91% sensitivity (95% CI = 0.86-0.94), 83% specificity (95% 
CI = 0.78-0.88), 83% positive predictive value (PPV) and 91% 
negative predictive value (NPV). 
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Validation of hCGres-predictive ability using the test data 
set. Individual values of hCGres were estimated in Test data set 
patients using the same model as defined above. Patient risk of 
resistance was predicted based on the > 20.44 IU1~^ hCGres 
threshold. The outcome regarding the predictive value of modelled 
hCGres in the Test patient data set was almost identical to that 
observed with the Model set patients: AUC = 0.87, (95% = 0.65-1) 
with 88% sensitivity (95% CI = 0.82-0.93), 86% specificity, (95% 
CI = 0.81-0.90), PPV = 79% and NPV = 92%, confirming the high 
reproducibility of the methodological approach. 

Comparison of hCGres against other predictors of MTX 
resistance in the combined multivariate analysis data set 
including 669 patients with required data 

Logistic regression. Using a binomial logistic regression model, 
including the WHO-FIGO score (HR= 1.17, 95% CI = 0.99-1.38), 
choriocarcinoma (yes vs. no, HR= 1.02, 95% CI = 0.36-2.87), high 
hCG titre at week 1 categorised by >2000IUl~^ cut-off value 
defined by Growdon et al (2009a; HR=1, 95% CI = 1-1), high 
hCG at week 7 categorised by >737IUl"^ threshold defined by 
Kerkmeijer et al (2009; HR= 1, 95% CI = 1-1.01) and hCGres, the 
only significant independent predictive factor of MTX resistance 
was hCGres (HR= 1.009, 95% CI = 1.007-1.011, P<0.0001). 

ROC curve analyses. The AUCs of the discrete ROC curves for 
other previously reported predictors of the MTX resistance in 
multivariate analysis data set were also assessed (Figure 4B). The 
highest AUC was with hCGres (> 20.44 lUl"^ AUC = 0.86, 95% 
CI = 0.83-0.89) compared with hCG value on week 1 
(>2000IU1"^ AUC = 0.64, 95% CI = 0.60-0.69), hCG values 
on week 7 (>500IUr^ AUC = 0.70, 95% CI = 0.65-0.74; 
> 520.24 lUl"^ AUC = 0.69, 95% CI = 0.65-0.74; >737IUl"^ 
AUC = 0.66, 95% CI = 0.62-0.70) and the WHO-FIGO score 
(AUC = 0.66, 95% CI = 0.62-0.70). 

Prediction of MTX failure. After a median follow-up period of 
4.89 years, 358 patients out of 800 patients (44.7%) experienced 
resistance and/or a relapse (i.e., MTX treatment failure). Relapses 
were observed in 19 out of 461 patients considered sensitive to 
MTX (4.1%) after a median of 277 days (IQR: 187-462). All 
patients with MTX resistance or relapse were treated with second- 
line treatment (Table 3). Three-year MTX failure-free survival was 
significantly higher in patients with hCGres values ^ 20.44, (88.4% 
vs 17%, P< 0.0001; Figure 4C). 

Data from 669 patients in the multivariate analysis data set was 
also analysed in a multivariate Cox model. The only significant 
independent predictors of MTX failure-free survival were as 
follows: high-modelled individual hCGres categorised by > 20.44 
threshold (HR= 14.51, 95% CI = 9.78-21.09); high hCG titre at 
week 1 categorised by >2000IU1~^ cut-off value defined by 
Growdon et al (2009a; HR = 0.62, 95% CI = 0.45-0.83) and high 
hCG at week 7 categorised by >737IU1~^ threshold defined by 
Kerkmeijer et al (2009; HR = 2.19, 95% CI = 1.65-2.90). 



discussion! 



Current trophoblastic disease management is heavily reliant on 
serum hCG monitoring (Seckl et al, 2010). As such, hCG and low- 
risk GTN provides an excellent disease model to demonstrate the 
potential utility of population kinetic modelling of disease-specific 
biomarkers for rationalisation of treatment decisions. 

Using the largest cohort of low-risk GTN patients reported 
in the literature to date, we have demonstrated that population 
kinetic modelling of patients' hCG measurements during the first 50 
days of MTX treatment provides a parameter, hCGres, which is a 
strong, independent and reproducible predictor for MTX resistance. 
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1 10 100 1000 10000 100000 0 10 20 30 

Individual predicted hCG values (lU 1"^) Time (days) 



Figure 2. Internal qualification of the model in Model data set patients (n = 418 patients). (A) Observed vs individual predicted hCG titres (blue 
circles), identity line (black line), smooth line of the cloud of points (red line). (B) Visual predictive checks of hCG concentrations. The observations 
are symbolised by the blue circles, and the median, 5% and 95% of observations are represented by red lines. Blue areas represent the 95% 
prediction interval for the median, 5% and 95%, computed fromlOOO simulated replicates. A full colour version of this figure is available at the 
British Journal of Cancer journal online. 




5 10 15 20 25 30 35 40 45 50 
Time (days) 



5 10 15 20 25 30 35 40 45 50 
Time (days) 




- Observed hCG 

Individual 
predicted hCG 



0 5 10 15 20 25 30 35 40 
Time (days) 



0 5 10 15 20 25 30 35 40 45 50 
Time (days) 



Figure 3. Examples of four typical patients of Model data set. The blue dots are the observed hCG values while the green line represents the individual 
predicted values. Resistances to MTX are specified. A full colour version of this figure is available at the British Journal of Cancer journal online. 



The continuous ROC AUC for discriminating subsequent MTX 
resistance in the Model data set was 0.94, which, according to Swets' 
(1988) criteria, describes an 'excellent' test. These findings are 
consistent with those of the first French study, and those of a recent 
retrospective analysis of American GOG- 174 phase III trial data (You 
et al, 2010; You et al, 2012). Moreover, the risk of MTX resistance 
could be predicted with better accuracy than using other previously 
reported predictors, including absolute hCG concentration thresholds. 
The kinetic approach was expected to have greater predictive power 
and show better reproducibility than using single hCG concentration 
thresholds as the effects of measurement intra- or interindividual 
inaccuracies (e.g., related to the use of different immunoassays) are 
reduced during fitting by the least squares method across multiple 
measurements (You et al, 2010; Bidart et al, 1999). 

However, the application and comparability of the cut-off values 
derived by the other studies is limited due to differences in (i) study 
design and categorisation of patients (Kerkmeijer et al (2009) 
study's derived cut-off value was for discrimination of 'mono- 
therapy resistance' rather than just MTX resistance alone), (ii) 
prevalence of the test positive condition and (iii) the use of 
different assays (van Trommel et al, 2006; Savage et al, 2008). 



The proportion of MTX-resistant patients in our low-risk 
GTN cohort (42%) was higher than the varying proportions 
reported in earlier clinical studies (ranging 18 to 40%; Powles et al, 
2007; Patel and Desai, 2010), but is in agreement with the largest 
studies which found 39-43% (Kerkmeijer et al, 2009; Sita-Lumsden 
et al, 2012). 

Most strikingly, when the same cut-off value derived for hCGres 
in the present study was applied to the original French study's data 
set of 154 patients, where serum hCG was measured weekly using 
various assays at different centres, an equivalent ROC assessment 
for MTX resistance was found with AUC = 0.85 (95% CI = 0.78- 
0.92), further confirming the high reproducibility of the metho- 
dological approach. 

The hCGres parameter also appeared to be a strong indepen- 
dent predictor of MTX failure, defined as MTX resistance or 
relapse, sometimes occurring several months after hCG normal- 
isation. Prediction of MTX failure could have greater clinical 
relevance than MTX resistance alone, as it may be argued that the 
better measure of a specific treatment's success is lifelong cure 
from the disease (Powles et al, 2007). Patients with an hCGres 
> 20.44 lU 1 ~ ^ after three fully dosed treatment cycles of MTX had 
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ROC curves for prediction of MTX resistance MTX failure free-survival 




Week 1 hCG > 2000 lU r'' 

Modelled hCGres > 20.44 lU 

WHO-FIGO score 

Reference line 

Figure 4. Predictive value of modelled kinetic parameters. (A) Predictions of MTX resistance offered by two modelled kinetic parameters K 
(AUG = 0.69, 95% CI = 0.64-0.74) and hCGres (AUG = 0.94, 95% CI = 0.91-0.96). (B) Predictive values of hCGres ( > 20.44 lU I " ^ : AUG = 0.86, 95% 
GI=0.83-0.89) against previously reported predictors: hCG values on week 7 (>737IUI"^: AUG = 0.66; 95% GI=0.62-0.70; >520.24 lU I " ^: 
AUG = 0.69, 95% 0=0.65-0.74; >500IU|-^: AUG = 0.70, 95% 0 = 0.65-0.74) and the WHO-FIGO score (AUG = 0.66, 95% 0 = 0.62-0.70). 
(C) Predictive value of modelled hCGres categorised by 20.44 cut-off value regarding MTX failure-free survival. Three-year MTX failure-free 
survival = 89.7% vs21.8%, hCGres < 0.0001. 



Table 3. Summary treatment and relapse details for patient cohort 







Response to MTX 




Sensitive 
(n = 461) 


Resistant 
(n = 339) 


Median MTX cycles received 
(range) 


7 (3-14) 


4 (3-11) 


Relapse 


No 

Yes 


442 (95.9%) 
19 (4.1%) 


ND 
ND 


Second-line treatment 


ACT-D 

EMA-CO 

Hysterectomy 


0 (0%) 
17 (89.5%) 
2 (10.5%) 


164 (48.4%) 
192 (51.6%) 
0 (0%) 


Median time to second-line 
treatment (range) 


1 1 .36 months 
(5.7-33.3) 


1 .8 months 
(0.9-5.2) 


Abbreviations: ACT-D = actinomycin D; EMA-CO = etoposide, MTX, ACT-D, cyclopho- 
sphamide and vincristine regimen; MTX = methotrexate. 



a more than 80% chance of MTX- resistant disease and a 92% 
3 -year MTX treatment failure risk. 

Our fmdings suggest a possible rationale to consider changing to 
second-line treatment immediately following prediction of MTX 
resistance using the kinetic method. In the present data set, if patients 
with hCGres values > 20.44 IU1~^ were changed to second-line 
treatment after only three fully dosed cycles of MTX treatment, a total 
of 478 MTX cycles administered to 226 patients, who were later 
determined to be MTX resistant and treated with ACT-D {n = 109) 
or EMA-CO (n = 117) regimens, would have been saved (83% true- 
positive rate). Such a strategy would reduce the duration of MTX 
therapy by 2 weeks in 50% patients (and by 4 weeks in 25% patients) 
and limit unnecessary use of MTX. However, this approach would 
also have led to unnecessary treatment change in 66 patients who, 
despite their high-modelled hCGres values, were cured with MTX 
therapy (17% false-positive rate). The strategy appears a fair one as 



only 9 out of 66 (14%, or 1%, of the entire cohort) of these patients, 
whose hCG levels where ^300IU1~^ after three courses of 
MTX, would, according to current clinical practice, have been 
switched to the higher toxicity EMA-CO regimen (Seckl et al, 2010), 
with minimal probability of mortality as recently shown (Alifrangis 
et al, 2013). The remaining 57 patients would have been changed 
to the lower toxicity ACT-D regimen, expected to be highly curative 
in such patients, as shown by the GOG- 174 trial outcomes 
(Osborne et al, 2011). Although unnecessary treatment with 
EMA-CO represents a greater clinical risk than prolonged ineffective 
treatment with MTX, where the agreeable balance between 
them exists is a matter for further research. The authors acknowl- 
edge that implementation of a strategy requiring weekly measure- 
ments of serum hCG titres may not be possible in all parts of the 
world. 

In summary, we have shown that of the previously reported efforts 
to utilise weekly hCG measurements to predict MTX resistance in 
low-risk GTN, the population kinetic modelling approach had the 
highest predictive power and was the most reproducible. We have 
shown that by using individual calculated hCGres values, it was 
possible to accurately identify after only three fully dosed MTX cycles, 
^ 80% of patients who will develop MTX resistance. Although model 
building is a complex process, our results may have practical clinical 
application. An online programme integrating the model could easily 
calculate individual modelled kinetic parameters, such as hCGres, 
from a patient's hCG results. We are currently designing a prospective 
multicentre European clinical trial with the European Organization 
for Treatment of Trophoblastic Disease, which would employ such a 
programme to further validate the clinical utility of hCGres for 
predicting MTX resistance, particularly with regard to reproducibility 
of predictive performance, inter-assay hCG cut-off values and the 
minimum interval requirement for testing. 
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